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Abstract 

We study numerically the SU(2) Kazakov-Migdal model of 'induced QCD'. In con- 
trast to our earlier work on the subject we have chosen here not to integrate out the 
gauge fields but to keep them in the Monte Carlo simulation. This allows us to measure 
observables associated with the gauge fields and thereby address the problem of the lo- 
cal Z2 symmetry present in the model. We confirm our previous result that the model 
has a line of first order phase transitions terminating in a critical point. The adjoint 
plaquette has a clear discontinuity across the phase transition, whereas the plaquette 
in the fundamental representation is always zero in accordance with Elitzur's theorem. 
The density of small Z2 monopoles shows very little variation and is always large. We 
also find that the model has extra local U(l) symmetries which do not exist in the case 
of the standard adjoint theory. As a result, we are able to show that two of the angles 
parameterizing the gauge field completely decouple from the theory and the continuum 
limit defined around the critical point can therefore not be 'QCD'. 



1 Permanent address 



Introduction 



Recently Kazakov and Migdal (KM) JT| presented an entirely new approach to the problem 
of finding a solution to QCD in the large N limit. In a nutshell, their philosophy is as 
follows: Suppose there is substructure to the quarks and gluons which becomes visible at 
energies large compared to some scale A. Then, regardless of the detailed nature of the 
constituents one should, upon integrating them out end up with QCD at the scale A as the 
only renormalizable, asymptotically free theory in four dimensions. KM propose to use a 
scalar field in the adjoint representation as a gluon constituent. The scalar mass acts as an 
effective ultraviolet cutoff for the theory. KM consider the following model, defined on a 
d-dimensional hypercubic lattice with action 



s = nj: 



(1) 



where <&(x) is a traceless scalar field in the adjoint representation of SU(N), covariantly 
coupled to the gauge fields U^{x) defined on the links of the lattice. At distances on the 
order of the lattice spacing this theory clearly looks very different from ordinary QCD defined 
on the lattice. But perhaps, KM argue, if we can find a critical point of the lattice theory 
at which the continuum limit can be taken and if furthermore we can manage to keep the 
scalar field heavy in that limit then QCD might be 'induced'. 

Let us see how the above scenario could possibly work in the limit of smooth weak gauge 
fields. Integrating out the scalar fields in the partition function will give rise to an induced 
action for the gauge fields which because of the gauge invariance of Eq. (1), will be given as 
an (infinite) sum over Wilson loops on the lattice. Moreover, because the scalar field is in 
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the adjoint representation, the Wilson loops will also be. In the case of V(<j>) = ^ptr^ 2 , the 
integral is Gaussian and one obtains the induced action ffl: 



2Kr)-i| tr[ /(r)| 
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l[T]m, 



where l(T) is the length of the path T. To leading order this model has a critical point 
at mQ = 2d as can be seen easily from Eq. (1) where the lowest order action becomes an 
action for a free scalar field. Expanding in fluctuations a F 2 V term is generated and using 
the definition m 2 = m 2 , — m 2 KM [|IJ obtained the correspondence 

This result arises simply from the scalar loop. In the scaling region this model is supposed 
to be QCD with m acting as an ultraviolet cutoff. Denoting by M g a physical glueball mass 
this immediately gives the usual relation between a physical quantity , the cutoff and the 
bare coupling 

m 2 a 2 487r 2 

~ TlNgl ' [ ) 
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Combining this expression with that of Eq. ([J) one gets a power law scaling 

M] = (m 2 r , (5) 

with 7 = 23/22! 

If we assume that the above arguments are in fact correct then the possible payoff of 
the particular choice of inducing field is indeed enormous: Based on the work of KM and 
of Migdal it looks as though in the limit of large N the model is analytically tractable 
and that the degrees of freedom associated with the eigenvalues of $ could actually be the 
long sought "master field" of large iV QCD. The model is easier dealt with analytically and 
numerically if one first integrates out the gauge fields. Migdal || has found scaling solutions 
to the equation for the eigenvalue distribution in the model and obtained 7 = 1.365.. in four 
dimensions. More recently he even succeeded in obtaining the spectrum of small fluctuations 
around these solutions^. 

To summarize this section: If the model defined by Eq. (1) is to induce QCD correctly 
we expect the system can be tuned to a critical point and power-law scaling of physical 
quantities with non-trivial (non mean field) critical exponents. 

The SU(2) Phase Diagram 

The road to proving that the model in Eq. (1) induces QCD is a long one. First one must 
demonstrate that the model for arbitrary N has indeed a critical point at which a continuum 
limit can be constructed. This is not a trivial point: Lattice Higgs models typically have first 
order phase transitions and critical behavior is more an exception than the rule. The KM 
model is actually not quite unknown to lattice theorists - in the case of SU(2) it is the lattice 
version of the Georgi-Glashow Q model in the limit of infinite gauge coupling. This model 
has been studied before but the phase structure in the case of infinite gauge coupling 
was until recently unknown. Once the existence of a critical point has been established one 
should (ideally) calculate the spectrum and check whether (5) is satisfied. Finally, in order 
for an exact solution of the theory at large N to be relevant for physics at N = 3 one must 
show that there are no unexpected phase transitions in N. As to the last point, the K-M 
model has recently been solved in one dimension at arbitrary iV and it was found that the 
TV dependence is smooth ||. However in 1-D the model is trivial unless one compactifies 
the line to a circle and even then it is effectively a one point model since one can always go 
to a gauge in which all the links but one are equal to unity. This seems like a rather severe 
over-simplificat ion . 

In the following we will study N = 2 case which is the simplest non-trivial version of the 
model. Our objective is to determine whether it is possible to construct a continuum limit 
for the theory and in what sense it could be equivalent to the 577(2) gauge theory ('QCD') 
in that limit. First we must search for a critical point. In our previous work || we did just 
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that by integrating out the gauge fields following KM's suggestion. One obtains 



/ [d<f>] exp J V \lnct) 2 - 2V{(j>)] + V In 
J <t»o I x J <xy> 



sinh(4(f)(x)4>(y)) 
<f>(x)(f>(y) 



(6) 



where < xy > denotes the sum over nearest neighbors of x. The field <f>{x) in Eq. 
is the (positive) eigenvalue of the matrix $. In this form the problem is to determine 
whether the above scalar theory has a nontrivial continuum limit. It is well known that 
the continuum limit of a scalar model with polynomial interaction terms is trivial with 
vanishing renormalized coupling J7J. In principle the continuum limit of Eq. (|5|) could avoid 
the triviality problem due to the presence of the nonpolynomial interaction terms. But note 
that these terms appeared as a consequence of the fact that we could integrate out the phase 
of the scalar by doing the U— integral exactly. If we could integrate out the angular variables 
in an O(N) model exactly the resulting effective action for the radial variables would also 
be nonpolynomial. The phase diagram that was found in Ref. H, corresponding to a scalar 
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potential of the form V(<f>) = '^f-trcf) 2 + |tr0 4 , is shown in Fig. 1. There is a first order 
phase transition at small A which ends in a critical point. In a simple mean field theory 
the critical value is A = 2.57, Wq = 4.52. From the Monte Carlo the precise location of the 
endpoint is very hard to determine. But it's existence is fortunate: It allows us to construct 
a continuum limit. Interestingly though, in the A — > limit no continuum limit exists. This 
is consistent with the observations by Gross || who constructed an exact large N solution 
in the case of a quadratic potential. 



The Nature of the Phase Transition 



In this section we will present the results of a Monte Carlo study of the theory as described 
by the action in Eq. (1). The gauge invariant content of the $- field lies in it's eigenvalues, 
so we again diagonalize $ and write 



Z = J [dU] J [d<f>] exp | Y, H 2 - 2V(0)] + 4 J2 <Kx)<K* + A*)(2|*»| 



-1^ 



(7) 



Here we have made use of the parameterization U = ( _ x _ J . Keeping the gauge fields 

in place will allow us to study the nature of the phase transition that was discussed in the 
previous section. 

Note that the theory has in addition to the local 577(2) gauge invariance additional local 
U(l) symmetries. This very important fact was pointed out by Gross ||. The center symme- 
try Un(x) — > Z^(x)U^(x) discussed by Kogan, Semenoff and Weiss [[[(]], which is responsible 
for the fact that only adjoint loops appear in Eq. (2), is a subset of this larger symme- 
try. We can make the larger symmetry very explicit by using ordinary gauge invariance to 
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diagonalize the $ field as in (7) to obtain the gauge field action (in unitary gauge') 



S G = 2^ (t>{x)(j){x + ^trlnU^x^U^x)] (8) 

X,fl=l 

The remaining symmetry now after the diagonalization is (V, W are diagonal matrices G 
577(2)) Ufj,(x) — > V At (x)f/ At (a:)W / At ^(x). Since one of them is the diagonal part of the usual 
SU(2) gauge transformation, we have (1x2-1 extra U(l) symmetries |. Now use the 
following (Euler) parameterization of the links 

v e-^ 1 ) { -b VT=¥ ) { e- ie > ) ' iyj 
In terms of these parameters we obtain 

d 

S G = AJ2 Hx)<p{x + fJ.)(2bf(x) - 1). (10) 

x,n=l 

Hence the reason for why (7) depends on only one real parameter associated with the links 
is apparent: When the $ field is diagonal the additional £7(1) symmetries can be used to 
eliminate two (diagonal) degrees of freedom from each link. Eq. (10) does not depend on 
the 9 angles in Eq. (9)|. Note that our counting differs by a factor of 2 from that of Gross, 
which apparently was based on counting degrees of freedom in the naive continuum limit. 
It is indeed true that in the naive continuum limit where = 1 + iA^ only the diagonal 
component A% drops out of the action. However as we can see from Eq. (10) the action 
depends only on one component of the gauge field. This is important: The naive continuum 
limit is in fact too naive in the KM model. The correct limit is obtained by writing the b 
dependent matrix in Eq. (8) as e zAr2 and expanding in A. The resulting action is (10) with b 
replaced by A. The significance of this observation lies in the fact that in the weak coupling 
(i.e. large 0) limit only the 6 M degree of freedom becomes smooth, the two angles 81,62 in 
Eq. (10) are always randomly distributed. The importance of the correct continuum limit 



was recently pointed out in a new publication of Kogan et al. [O . 



In this paper we have monitored in addition to the <fi field also the expectation value of 
the plaquette in the adjoint representation 

P* = \(\trV{P)\ 2 -l). (11) 

Due to the local Z2 symmetry of the model Wilson loops in the fundamental representation 
vanish according to Elitzur's theorem [|I2 . Based on the behavior of this quantity in addition 



to that of the 'density of Z 2 monopoles' 



p=<l- J] sign(trU(P)) > (12) 
Pec 3 



2 It is hard to see this symmetry in the induced action, eq.(2). 

3 In SU(3) we also found a parameterization in which U — D1VD2 (D\, D2 are diagonal matrices £ SU (3)) 
so that the action in this case only depends on 4 parameters. We suspect that this generalizes to SU(N): 
the action only depends on (N — l) 2 parameters. 
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we will later on speculate on the significance of the Z 2 degrees of freedom in the KM model. 
Eq. (9) and (10) imply that we can derive a simple expression for P a( ij- The Haar measure 
for the parameterization (9) is 

dU = 2bdbd6 1 d8 2 . (13) 

Using the fact that b enters the action only quadratically we can completely integrate out 
the gauge field to obtain a formula for the adjoint plaquette in terms of <fi. Denoting by 

C= ( 1 ~ a , a ), (14) 
V a 1 — a J 

where 

a = - , (15) 

a matrix on the linksQ which depends on at the beginning and the end of the link through 
uj = (p(x)(p(x + n) we obtain 

< P adj >=\< MC&CsCt) - 1) > . (16) 

The average is over the <fi field. In mean field theory, where we set <fi = const and determine 
the mean field from the minimum of the resulting potential one obtains the simple result 

<P adj >= l -(l-2a MF )\ (17) 

Later this will be compared to the result of our simulation. Note one important consequence 
of this result: The normalization of P a dj is such that in the naive continuum limit P a dj = 1. 
This is obviously wrong though. From Eq.(14) we see that in the weak coupling limit [] 
{uj — > 00) P ad j — > ~. Again this agrees with the observations of Ref. [pj] . It is also easy to 
see that the angular integrations imply that p — 1. 

We have simulated the theory described by the action in Eq. (|7|) on both 4 4 and 8 4 
lattices using a simple Metropolis algorithm on both the links and u where <ft = e u . It turns 
out that for the local quantities we measure there is hardly any difference in the results on 
the two lattices. Of course this is consistent with the fact that the transition is first order. 
As far as the the expectation value is concerned we reproduce the results of Ref. [§] where 
the gauge fields were integrated out exactly. This is a good check on the program. In Fig. 2 
and 3 we show a thermal cycle for and the adjoint plaquette, < ^(\trU(P)\ 2 — 1). There 
is a clear discontinuity in this quantity just as one would expect based on simple mean field 
theory. Clearly the agreement between Monte Carlo and mean field theory is excellent. The 
adjoint plaquette close to the critical point is shown in Fig. 4. The magnitude of the gap 
(if there is any) has decreased considerably. Also note that the adjoint plaquette rises only 
very slowly. At small m\ it is only 0(0.14) which is far away from the weak coupling limit 
Padj = §• Hence relatively close to the end point of the first order phase transition the gauge 
field is far from being perturbative! 



4 In general {U^) lj U k i -> 5u8j k (C)ij . 

5 In the strong coupling limit (uj — ► 0) P a dj — * 0. 
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We have also measured the fundamental plaquette in various places in the phase diagram 
both above and below the phase transition line in Fig. 1. It is always zero. An interesting 
quantity is the local plaquette which is histogramed in Figures 5 and 6. Note how in the 
weak coupling (small wIq) phase peaks develop at \trU = 1, —1. One might think at first 



that this has something to do with some kind of l Z 2 Higgs' transition advocated in Ref. flip 
However the explanation is much simpler than that: At weak coupling, when 6^0, the 
local plaquette is roughly given by 

^trU(P) « cos(6 P ) (18) 

where 8p = (81 + 82) p- Averaging over dp with uniform distribution gives a vanishing 
plaquette but the (normalized) distribution will be approximately 

W(-trU(P)) « — ; 1 =. (19) 

We have also measured p. This quantity is always large, p ~ 0.96 showing hardly any 
variation as one moves around in the phase diagram. As we mentioned before based on the 
randomness of the two angles the expected value of this quantity is unity. In the next section 
we will discuss the significance of this observation. 



Discussion and Conclusion 



In light of Eq. (2) one is tempted to view the KM model as a theory with a generalized 
action in the adjoint representation. The novel feature here is that there is actually an 
infinite sum over arbitrary shaped loops. In the naive continuum limit this action seems to 
be as good as any other in that it produces a term proportional to F^ 2 , with a non Abelian 
field strength F (the Abelian component drops out). But as we saw in the last section the 
naive continuum limit is incorrect in this model and we must be careful here. 

What is the mechanism that drives the phase transition? At first look the KM model 
reminds us of the studies of the 50(3) invariant theory 10 years ago. It was found that the 
mixed action theories have first order phase transitions in the (ftp, (3 a) plane, including the 
point f3p = [Ol]. An elegant explanation of why one expects a first order phase transition 



in the 50(3) theory was given by Halliday and Schwimmer [14|: The theory contains zero 
energy monopoles (in contrast to the SU(2) case where the presents of monopoles costs 
energy) which condense at some critical value of the coupling. In the monopole free phase 
the adjoint and fundamental theory approach a common continuum limit. Halliday and 
Schwimmer base their considerations on the following 50(3) invariant theory: 

5 = (3j2 a ( p ) trU ( p ) ( 20 ) 
p 

where the sum is over all plaquettes on the lattice. The variable a takes on values in Z 2 
and lives on plaquettes. Under Z 2 it transforms as o — > (HieP z i)°~- A string of plaquettes 
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carrying a nontrivial a can be thought of as the Dirac string of a monopole. The location 
of the monopole can be identified with an end of the string defined as a point on the lattice 
where 

p= n ( 2i ) 

Pec 3 

is negative. Here C3 is a three dimensional cube on the lattice. Halliday and Schwimmer 
argue that in the phase in which the monopoles are absent we can write a ~ ILgp z i m which 
case the z\ can be absorbed into the links. In particular for a 'tiled' Wilson loop we have 

< n^(o n °( p ) >u~< n >u> m) 

lec PeA(c) lec 

with U'(l) = ziU(l). The above equation tells us that 'tiled' (i.e. Z 2 invariant) Wilson loops 
in the adjoint theory described by Eq. ( p2j) will behave like ordinary Wilson loops in the 
theory defined using the fundamental representation. Note, that this does not mean that 
the fundamental Wilson loop is ever non zero in the adjoint theory. This is not allowed 
as a consequence of Elitzur's theorem [El. For this reason we do not quite understand the 



analysis in Ref. |10[ where it is claimed that there can be a 'Higgs' phase with a nonvanishing 
fundamental plaquette. Monopole condensation can be checked by computing the monopole 
density < p > which is invariant and always non-zero. Monte Carlo simulations [[ljj show 



a sharp rise of this quantity across the transition making the case for the above scenario of 
the phase transition quite convincing. 

Does this scenario apply to the SU(2) KM-model? We do not think so. We have seen 
that the quantity p is always large and shows absolutely no variation as we go through the 
phase transition. Hence Z 2 monopoles are always condensed in our system! This is of course 
due to the fact that as far as the 9 degrees of freedom are concerned the theory is always in 
the strong coupling phase. We therefore conclude that the first order phase transition is not 
a Z 2 driven transition. If one wanted to understand the phase transition entirely based on 
the formulation in terms of adjoint loops one must keep all the loops and probably introduce 
'large Z2 monopoles' for all of them. These 'monopoles' will interact and the picture will 
become a hopeless mess. Of course somehow this 'covariant' view (as opposed to 'unitary' 
one based on diagonalizing $) of the phase structure must be equivalent to the much simpler 
description of the phase transition in terms of Eq. (10). 

In the description of Eq. ( |IUD the Z 2 degrees of freedom have completely decoupled and 
are just part of the randomly fluctuating 9 variables. The transition is described very well 
by simple mean field theory which 'explains' the transition in the same way we explain for 
example the liquid-gas transition. We have seen that the one relevant degree of freedom of 
the gauge field becomes perturbative at small A (in the weak coupling phase). Close to the 
end point of the transition line the adjoint plaquette is still rather small. This is of course 
due to the fact that (f) is small there. Assuming for the moment that by carefully tuning 
the parameters of the scalar potential we can construct a continuum theory, what will it be 
like? From what we have seen it certainly does not look like the £77(2) gauge theory. Since 
two internal degrees of freedom have disappeared per component of the gauge field the most 
likely scenario is that we will just get a (trivial?) theory involving a massive photon (with 
two transverse and a longitudinal components) and a scalar. 
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For large N things could be different: There the theory retains it's non Abelian nature 
even after all symmetries have been used to eliminate degrees of freedom. Even in this case 
the question remains in what sense the theory becomes equivalent to QCD at the critical 
point. One would for example like to establish that there is a string tension in the continuum 
limit. Fundamental Wilson loops always vanish in the KM-model whereas adjoint loops have 
a perimeter law due to screening. Maybe the answer are the 'tiled loops' (Wilson loops whose 
minimal area is tiled with plaquettes) introduced by Kogan et al. 0] although they have 



the strange property of not vanishing at infinite coupling. Migdal in the meantime |Tb| has 
proposed to discard the simplest version of induced QCD and instead add in addition to 
the adjoint scalar also rif (1 << rif << N) fundamental fermions to the action. In this 
approach the fundamental Wilson loop need not vanish and may contain in addition to the 
perimeter piece also an area law contribution. We have also just become aware of a new 



paper by Kogan et al. [17 where it is argued in mean field theory (577(2)) that critical 



behavior could be recovered in the case A = by changing the measure for the scalar. While 
this might be true, based on what we have found the resulting theory will not be 'QCD'. 
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Figure Caption 



Figure 1: The phase diagram. The solid line is the mean field theory prediction and 
the simulation results are indicated by diamonds. The end point of the first order phase 
transition line is marked by * for MFT and square for Monte Carlo estimate. At the end 
point the theory becomes critical. 

Figure 2: Comparison of the mean field theory and Monte Carlo simulation result for 0-field 
vacuum expectation value. The solid lines are the MFT predictions and Monte Carlo data 
points are indicated by *. The errors are smaller than the size of the symbols. 

Figure 3: Comparison of the mean field theory and Monte Carlo simulation result for the 
adjoint plaquette. The solid lines are the MFT predictions and Monte Carlo data points are 
indicated by *. 

Figure 4: Monte Carlo simulation result for the adjoint plaquette. The value of A is close 
to the critical value. 

Figure 5: Histogram of the local plaquette in the weak coupling phase. 
Figure 6: Histogram of the local plaquette in the strong coupling phase. 
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